Enhancement of magnetic anisotropy barrier in long range interacting spin systems 
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Magnetic materials are usually characterized by anisotropy energy barriers which dictate the time 
scale of the magnetization decay and consequently the magnetic stability of the sample. Here we 
consider magnetization decay in a class of spin systems in contact with a heat bath, with long range 
interaction and on site anisotropy. We show that anisotropy energy barriers can be determined from 
ergodicity breaking energies of the corresponding isolated system, independently of the mechanism 
of magnetic decay (coherent rotation or nucleation). In particular, we show that, in presence of 
long range interaction the anisotropy energy barrier grows as the particle volume, V, for sufficiently 
large volume, and it can grow, for finite sized systems, as y 2 - a / d ^ w h e re a < d is the range of 
interaction and d is the embedding dimension. As a consequence there is a relevant enhancement 
of the anisotropy energy barrier compared with the short range case, where the anisotropy energy 
barrier grows at most as the particle cross sectional area for large or elongated particles. 

PACS numbers: 05.20.-y,05.10.-a, 75.10.Hk, 75.60.Jk 



I. INTRODUCTION 

A truly comprehensive understanding of magnetism at 
the nanoscale is still lacking. From a theoretical point of 
view the problem of magnetization decay in nanosystems 
is difficult to treat: nanoscopic systems are too big to 
be solved by brute force calculation and too small to 
be tackled by the tools of statistical mechanics at the 
equilibrium. Indeed, the problem of magnetization decay 
is a typical example of out of equilibrium phenomenon, 
which is the decay out of a metastable state. 

On the other hand, nanomagnetism has also impor- 
tant consequences in the technology of memory and in- 
formation processing devices. The quest for improv- 
ing magneto-storage density calls for the realization of 
smaller and smaller magnetic units. Significant improve- 
ments in experimental techniques allowed investigations 
of magnetic properties in nanoparticles and nanowirespQ. 
In particular, recently, there has been great interest in 
Single Chain Magnets (SCM) gHS], which are possible 
candidates for nanoscopic memory units. Nanoscopic 
systems can also show ferromagnetic behavior at low but 
finite temperature, even if a ferromagnetic phase transi- 
tion is theoretically forbidden [7], due to large magnetic 
decay times. 

The modelization of the magnetization decay, through 
an on-site anisotropic barrier is quite typical in litera- 
ture, even if, mainly short-range interactions have been 
considered. Nevertheless, in many realistic situations, 
one needs to go beyond nearest neighbor coupling, tak- 
ing into account the long range nature of the interaction 
defined by a two-body spin interaction constant decay- 
ing at large distance with a power law exponent a not 
larger than the embedding spatial dimension d[8j. It is 
the case, for instance, of the dipolar interaction in 3-d 
systems, or of the so-called RKKY (Ruderman-Kittel- 
Kasuya-Yosida) interaction, which decays as r~ d , where 



r is the distance between spins, and d is the dimension of 
the lattice system. In particular, the latter might be re- 
sponsible for the ferromagnetic behavior of Diluted Mag- 
netic Semiconductors (DMS) [5] and Diluted Magnetic 
Oxides (DMO) [10 , promising materials for the realiza- 
tion of spintronics devices. 

One of the first attempts to understand magnetic de- 
cay times in nanoparticles is due to Nccl[ll and Brown 
|12j , who considered that all the spins in a magnetic par- 
ticle move coherently, so that they can be considered as 
a single spin and described magnetization decay as due 
to thermal activation over a single energy barrier. In 
Brown's theory this time, r, is shown to follow an Arrhe- 
nius Law : 



where /3 = l/fc^T is the inverse temperature and AE oc 
V is the anisotropic energy barrier proportional to the 
particle volume V . A step forward Brown's theory has 
been realized by Braun |13j . In his theoretical approach 
a sufficiently elongated system of short range interact- 
ing spins with an on-site anisotropy barrier, have been 
shown to reverse their magnetic moments (thus produc- 
ing an average magnetization decay) through a process 
called nucleation, energetically convenient with respect 
to coherent rotation. In this mechanism, accomplished 
by the formation of a soliton-antisoliton domain wall, 
the magnetic anisotropic energy to be overcome turns 
out to be proportional to the cross sectional area of the 
particle, AE cx A and not to its volume V. Studies 
of different mechanisms of magnetic decay have been the 
objective of intense investigation|14| until recently, where 
also 3-d spherical samples with short range interactions 
and on-site anisotropy are shown to produce nucleation 
for sufficiently large radius [15]. Thus, for short range in- 
teraction, Brown's theory, and a consequent Arrhenius 
Law with an exponent proportional to the volume V of 
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the particle is valid only for very small particles, while in 
general, for large or elongated particles, the exponent is 
given by the cross sectional area of the particle. A smaller 
exponent means smaller decay times for the same tem- 
perature. The size and shape dependence of the magnetic 
anisotropy barrier, and consequently of the decay times, 
have been also confirmed experimentally in Ref. [16 . 

The consequences of long range interaction on mag- 
netic decay is much less investigated. Long range in- 
teraction can affect the decay out of a metastable state 
in a significant way; in the seminal paper [IT] it was 
shown that the decay time out of a metastable state in 
a toy model with infinite range interaction is given by 
the Arrhenius Law with an exponent proportional to the 
squared volume of the particle. 

The main goal of this paper is to analyze magnetic 
decay beyond nearest-neighbor interaction, focusing on 
realistic long range interacting systems. In order to esti- 
mate the dependence of magnetic decay times on temper- 
ature, we propose a different point of view, which turns 
out to be independent of the decay mechanism and re- 
lated to the recently found Topological Non-connectivity 
Threshold (TNT) in anisotropic spin systems [HI [19] . 

II. TOPOLOGICAL NON-CONNECTIVITY 
THRESHOLD AND MAGNETIZATION DECAY 
TIME 

In this Section, following Ref. jT5], we briefly review 
the Topological Non- Connectivity Threshold. 

Let us consider a generic anisotropic spin system, with 
an easy axis of magnetization, (the direction n easy of the 
magnetization in the ground state), with a microcanoni- 
cal energy 

H(Si, . . . , Sn) = E. 

Let us also set 

m = (1/AQ }^ Sk ■ Ueasy, (2) 
k 

as the magnetization along the easy axis. Note that in 
our paper it will be n easy = z. 

It was proven [TH] that below a suitable threshold, E tnt , 
given by the minimal energy attainable under the con- 
straint of zero magnetization m along the easy axis: 

E tnt = Mm(H(...Si...) | m = 0), (3) 

the constant energy surface is disconnected in two por- 
tions, characterized by a different sign of the magnetiza- 
tion. From the dynamical point of view one has a case 
of ergodicity breaking: a trajectory at fixed energy can- 
not change the sign of magnetization since it is confined 
forever in one region of the phase space. 

It was also demonstrated that in case of long range 
interaction among the spins ffiT), the disconnected energy 



portion determined by the TNT, remains finite, when the 
number of particles becomes infinite. 

While for isolated systems, the magnetization cannot 
reverse its sign if its microcanonical energy E is below the 
energy threshold E tnt , when the system is put in contact 
with a heat bath this may happen for any (even extremely 
low) temperature. Question arise whether the TNT dic- 
tates in some way the magnetization decay. This is ex- 
actly what we have found (see Sections below), namely 
AE tn t = Etnt — E m i n , represents an effective energy bar- 
rier for magnetic decay and the decay time depends ex- 
ponentially on such energy barrier: 



where, r is a factor, that may depend on temperature 
too, E mm is the ground state energy and j3 = 1/fcsT. 
While this was confirmed in simple toy models with all- 
to-all interaction among the spins [21] , here our aim is to 
generalize these results to realistic spin systems. 

It is possible to give a heuristic justification of Eq. Q. 
Magnetic decay occurs through fluctuations of the mag- 
netization around its equilibrium value. The probability 
a magnetization fluctuation is determined by the free en- 
ergy barrier, AF — AE—ksTAS, through the Arrhenius 
factor, e~P AF . Since the entropic barrier, AS, is usually 
negligible at low temperature, the accessible spin config- 
urations can be determined minimizing the energy only. 
In order to reverse its sign, the value of the magnetization 
has to go, say, from a magnetization m = 1 to m = 0. 
Since for m = 1 the system is in its minimal energy, it is 
clear that AE tn t represents the minimal energy barrier 
found by the system while reversing its magnetization. 

III. THE MODEL : ISOTROPIC a RANGED 
PLUS ON-SITE INTERACTION 

Let us now focus on spin systems with isotropic long- 
range exchange interaction and on-site anisotropy, de- 
scribed by the following Hamiltonian: 

where, Si are the spin vectors with unit length, a de- 
termines the range of the interaction among the spins, 
J > is the exchange coupling and D > is the on-site 
energy anisotropy. The minimal energy for this class of 
spin systems is attained when all the spins are alligned 
along the z direction, which thus defines the easy axis of 
the magnetization. 

In the following we will mainly focus on the "critical" 
case a = d, since it is relevant for realistic applications 
(dipole interaction in 3-D and RKKY interactions); here 
critical is referred to the fact that interaction is defined 
as long-ranged for a < d, while it is short ranged for 
a > d 8 . Moreover, we choose Hamiltonian ^ for sake 



of generality since we would like to consider not a single 
specific interaction, but a class of interactions depending 
as some power of distance among spins. 

For this class of systems, the energy E tn t, see Eq. Q, 
can be computed numerically using a minimizing con- 
strained algorithm. We can also estimate analytically 
E tnt for a generic range of the interaction a. To this 
purpose let us consider two configurations with m = 0: 

• the first one with all spins aligned perpendicular to 
the easy axis. The energy difference of this configu- 
ration from the one having minimal energy is DN, 
which is the energy barrier due to the coherent ro- 
tation of all spins; 

• a configuration, labelled fj,, consisting ot two 
neighbors identical blocks with opposite magnetiza- 
tion along the easy axis. This configuration roughly 
corresponds to what is called nucleation configu- 
ration in literature. The energy difference of this 
configuration from the minimal energy in the case 
ot < d, is given by [20] : 

AE n = JC a , d N 2 -^ d , (6) 

where C Qj( j is a suitable constant, for instance, 
C* M ~ 4 In 2. 

Following similar reasoning as in Ref. [20], it can be 
shown that the energy of these two configurations is a 
good approximation of E tn t, so that we can write 
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AEtnt = E f , 



E„ 



Min(DN,AE n ), (7) 



Note that Eq. ([7]), valid whenever DN is not close to 
AE^, gives an estimation of the anisotropic energy bar- 
rier, which can be used in Eq. |4) to get magnetic decay 
times. Eq. ^ has been confirmed numerically in Fig. (JlJ, 
for the case a = d = 1, and we have found coherent rota- 
tion for D <C D cr = 4Jln2 and nucleation for D 3> D cr , 
where the D cr has been obtained comparing Eqs. ^ and 
[H]. In this fi gure we compute numerically the energy 
barrier AE tnt using a minimizing constrained algorithm 
(symbols) and we compare it with the approximation ([7]) 
for one dimensional chains with different number of spins. 
As one can see, the agreement is fairly good. Note that 
in case of nucleation (DN > AE^), the effective energy 
barrier AE tnt to be put in the Arrhenius Law should be 
increased by a term D, called single-spin flip. This extra 
term should be added only for continuous models, as the 
modified Heisenberg model we are considering, so that, 
t oc e D " [3] . This should be done since a single spin can 
continuously change from S = 1 to S = —1 only pass- 
ing the state S — 0. While the two extreme states have 
the same magnetic anisotropy (since it is proportional to 
DS 2 ) the intermediate state 5 = has a barrier higher 
by a factor D. This additional spin-flip docs not occur, 
for instance, in the Ising-like models, since there are no 
intermediate states between 5 = 1 and S = — 1. 

In the following we analyzed the magnetic decay time 
in the canonical ensemble, using a modified Monte Carlo 
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FIG. 1: Energy barrier as a function of the anisotropy barrier 
D for the a = d = 1. Different symbols indicate different 
number of spins in the linear array, as indicated in the legend. 
Dashed lines indicate the approximation Eq. Q. The vertical 
black line indicates D = D cr , and separates the region of 
coherent rotation (CR, left), from the region of nucleation 
(NU, right). 



simulation [TU [52]. As initial condition we choose all 
spins aligned along the easy axis, and from the exponen- 
tial decay of the average magnetization, (m(t)} oc e -t / T , 
we computed the magnetic decay time, r. 



IV. ENHANCED BARRIER 

In this Section, we show our numerical results for the 
magnetization decay. They will clearly give evidence of 
the enhancement of the magnetization decay time due to 
long range interaction. 



A. One dimensional case 

Here we focus on the one dimensional critical case. 
Our results for the magnetic decay time, are collected in 
Fig. [2| where we consider the case a — d — 1: magnetic 
decay times are shown in the upper panel for coherent ro- 
tation (D <C D cr ) and in the lower panel for nucleation 
(D D cr ). Samples with different number of particles 
experience different anisotropy energy barriers. As one 
can see the numerical results indicated by symbols well 
agree with the theoretical prediction given by Eq. |7|, 
shown as full lines. Our results for coherent rotation, 
Fig. [2] (upper panel), clearly indicates that Brown's the- 
ory of coherent rotation with an anisotropy energy barrier 
proportional to the particle volume (V oc N ) is still at 
work even for long range interacting systems. In Fig. [2] 
(lower panel), the case of nucleation is shown. Also in 
this case the anisotropy energy barrier is proportional to 
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FIG. 2: Magnetic decay times, r vs the rescaled inverse tem- 
perature, /3/Pstat for a = d = 1 and different numbers of spins 
N as indicated in the legend. Upper panel is characterized 
by coherent rotation with parameters J = 1 and D = 0.05 , 
while the lower by nucleation with J = 1/16 and D — 0.5. 
Different symbols refer to numerical results, while full lines 
in are the theoretical predictions exp(/3 AEtnt)- Specifically 
for coherent rotation, we have AE tn t = DN, while for nucle- 
ation, we have AEtnt = AE-f^ + D (see text for explanations). 
Vertical dashed lines refer to the inverse rescaled statistical 
temperature /3 = /3 sta t obtained by a standard mean field 
approach. 



the particle volume. Let us remark that this behaviour 
is different for short range interacting systems, where in 
the one dimensional case, the anisotropic energy barrier 
is independent from the size of the system [21 [B] . Note 
that for high temperature (small (3 values) symbols with 
different number of spins lye upon the same curve, that 
turns out to be almost independent of N. 



B. Multidimensional case 

Our results can be extended to higher dimensions. In 
particular we will show that for large enough particle, the 
anisotropic energy barrier is always proportional to the 
volume of the particle. Moreover we point out that for 
finite sized systems, the anisotropic energy barrier can 



grow even faster than the volume. 

From Eqs ^ and Q, we have that while for coherent 
rotation, the anisotropy energy is always proportional to 
the volume of the particle, for nucleation the anisotropy 
energy can grow even faster than the volume of the parti- 
cle, AE-fi oc N 2 ~ a / d . Since the anisotropy energy barrier 
is determined by the smallest of the two barriers, it fol- 
lows that for sufficiently large particles the energy barrier 
will always be proportional to the particle volume since 
V oc N. In particular, this means that even in the case 
of 3-d elongated particles, differently with what happens 
in the short range case, the energy barrier does not de- 
pend on the cross sectional area but on the whole particle 
volume. 

This theoretical prediction has been numerically con- 
firmed in Fig. [3j where we study the case of nucleation in 
thin parallelepipeds LxLx eL, in presence of a " critical" 
interaction a = d = 3. We chose the aspect ratio e = 10. 
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FIG. 3: Energy barrier AE as a function of the side L of the 
parallelepiped, obtained by an exponential fitting of the decay 
time vs the inverse temperature /3. Here is a = 3, J = 1/20 
and D — 0.5 (nucleation). 

The energy barrier AE extracted from the numerical 
fit of the Arrhenius Law r = r exp(/3AE) has been plot 
as a function of the side L of the parallelepiped. For sake 
of comparison the lines proportional to the cross section 
(L 2 ) and to the volume (L 3 ) has been drawn. As one can 
see without any doubt the barrier grows as the volume. 

From Eqs Q and Q, it also follows that, when a < d 
the anisotropy barrier can grow even faster than the vol- 
ume of the particle. This can happen up to a critical 
value of the particle volume, for which DN = Ai?^. 
Such an effect is a distinguished feature of long-range in- 
teraction and it is forbidden in the case of short range in- 
teractions. To confirme such theoretical prediction let us 
consider a square lattice (d = 2) characterized by an ex- 
change interaction decaying as 1/r {a = 1). Even in this 
case we can distinguish between the regime of coherent 
rotation (large anisotropy barriers D) from that of nucle- 
ation (small D values) . The results of our simulation have 
been shown in Fig. |2] In the lower panel the theoretical 
predictions have been indicated as full lines while symbols 
represents numerical results from Montecarlo simulation. 
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FIG. 4: Decay time r us the rescaled inverse temperature 
/3/pBtat for a square lattice (d = 2) with a = 1. for J = 1 and 
D = 0.025 (upper panel, coherent rotation); and for J = 1/40 
and D — 1/2 (lower panel, nucleation). Symbols refer to 
numerical data, while full lines are the analytical prediction 
Eq. |7f. Vertical dashed lines represent the rescaled inverse 
statistical temperature P//3 s tat ~ 1. Note that in case of co- 
herent rotation (upper panel), the anisotropy energy is pro- 
portional to DN, while in the case of nucleation (lower panel), 
the anisotropy energy is proportional to JN 3 ^ 2 . Lattice di- 
mensions have been indicated in the legend 



As one can see the agreement with the theory, which 
predicts an energy barrier AE tnt oc J7V 3 / 2 is fairly good. 
For sake of comparison, the regime of coherent rotation 
(Fig. |4j upper panel) is also shown, indicating that, in 
this case, the energy barrier is proportional to the vol- 
ume of the particle, namely AE tnt oc DN. 

Another important point to be discussed is in which 
temperature range the Arrhenius Law with the exponent 
given by Eq. ^ holds. In the cases of coherent rotation 
and nucleation, see also Ref. [12], we might expect that 
the Arrhenius Law is valid only when kgT <C AE tnt . 
Clearly this gives an upper bound for the temperature for 
which Eq. Q is valid. Moreover it should be T < T stat , 
where the latter is the temperature at which the barrier 
at m = in the free energy vanishes (and it coincides 
with the temperature at which a phase transition occurs 
in the thermodynamic limit). It is very interesting that, 
computing T sta t = ^/ksPstat by means of a standard 
mean field approach, one gets a very nice estimate of the 
validity range of the Arrhenius Law (see dashed vertical 



lines in Figs. [2]and[4|. Nevertheless we stress that this 
delicate point needs further analysis to be clarified. 



V. COHERENT ROTATION AND 
NUCLEATION 



Let us notice that the use of the words coherent ro- 
tation and nucleation, are not only ways to label two 
different mechanisms taken from the short range models, 
but correspond to the spin dynamics. To prove that we 
first consider in Fig. [5] two different situations taken from 
the upper (coherent rotation) and the lower panel (nu- 
cleation) of Fig. [4] In particular we consider one single 
Montecarlo trajectory and we plot, in Fig. [5] the value of 

M = \ jf Sili @i\ an( i m = jf Dili Sii as a function of 
the Metropolis time. It is clear that while in the case of 
coherent rotation (upper panel of Fig. [5]) M s» 1 and m 
change its sign, for nucleation (lower panel of Fig. [5| one 
has M ~ \m\. 




100000 



100000 



FIG. 5: Upper panel (coherent rotation) for a square lattice 
6 x 6, P = 5. Lower panel (nucleation) for a square lattice 
8x8 and j3 — 4.5. In both panels the red curve refers to M, 
while the black one to m. 



As an interesting consequences of the coherent motion 
of the spins during the magnetic reversal, we show here, 
focusing on the one dimensional case, that one can esti- 
mate the whole probability distribution of the magneti- 
zation in the canonical ensemble, Pr(m) at low temper- 
ature and around m = 0. 

Note that the knowledge of the probability distribution 
of the magnetization allows one to compute the magnetic 
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decay time: 



= /3AF _ 



Pr{ni max ) 

PT(m m in) ' 



(8) 



where Pri^max^in) stands for the maximum/minimum 
value of the probability distribution. 

Due to the coherent motion of all spins it is possible to 
express the energy of the system as a function of the total 
magnetization along the easy axis; from the knowledge 
of E(m) we can compute Prim) since: 

P T (m) dm = P T {E) dE = e- E(m)/kBT p(E) dE 

where p(E) is the density of states at fixed energy. Us- 
ing the fact that p(E)dE = p(m)dm, where p(m) is the 
density of states at fixed magnetization, we have: 



P T (m) dm = e -E(m)/k B T p ( m ) dm _ 



(9) 



Let us compute E(m) for coherent rotation and nu- 
cleation. When the decay of the magnetization oc- 
curs through coherent rotation we have SiSj — 1, and 



S* = (S z ) = m so that, from Eq. ^ we have: 
E(m) = -J/2j2 1 / r tj ~ DNm 2 . 



(10) 



Since all spins move coherently on a sphere, the den- 
sity of states at fixed magnetization along any axis is 
constant, p(m) = const. Moreover the first term in the 
r.h.s. of Eq. ( 10 1 is independent of m and we have (for 
coherent rotations): 



P T (m) oc e? DNm2 . 



(11) 



We verified numerically that Eq. (11) gives an excel- 



lent approximation of the probability distribution of the 
magnetization around m = in the coherent rotation 
regime, see Fig. 6] (upper panel). 

In case of nucleation the possible configurations which 
the system can visit can be obtained assuming that mag- 
netic decay occurs by first reversing one of the edge spins, 
then its nearest neighbor, and then the spin immediately 
after until all spins are reversed. If we have k spins on 
the left pointing upwards and N—k spins pointing down- 
wards, we can write the energy of this configuration as 

E(k,N) = E min (k) + E min (N-k) + W(k,N-k), (12) 

where W(k 1 N — k) is the interaction energy between the 
two neighbors blocks of k and N — k spins and 

k k 

E mm (k) = -J/2 ]T 1 Kj ~ Dk > ( 13 ) 

is the minimal energy for a system of k spins. The mini- 
mal energy of a system of N spins can be easily obtained 
from ( 12 ) by changing the sign of W(k, N — k), 



Taking into account that m = (2k — N) /N is the magne- 
tization of a system composed by k spins along the easy 
axis and N — k opposite to it, one gets, 

E(m) = 2E min (N+) + 2E mm {N-) - E min (N), (15) 

where N± = N(l — m)/2. In this case it is more difficult 
to estimate p(m). Assuming, also in this case p(m) = 
const., from the knowledge of E{m) (Eq. ( |15[ )), we obtain 
Py(m) cx e~P E (™^ which agrees very well (apart close to 
m = ±1) with numerical results in the nucleation regime, 
see Fig. [6] (lower panel) . 




E min (N)=E min {k)+E min (N-k)-W(k,N-k). (14) 



FIG. 6: Probability distribution of the magnetization Pr(m) 
vs the magnetization m for a = d = 1, N = 20 and differ- 
ent temperatures j3 — 1/fcsT. Upper panel is for J = 1 and 
D = 0.05 (coherent rotation regime) and P = 6 (upper red 
hystogram), /3 — 9 (lower black hystogram). The lower panel 
is for J = 1/16 and D — 0.5 (nucleation) and f3 = 8 (upper 
red hystogram), /3 = 11 (lower black hystogram). Hystograms 
represent numerical data, while smooth curves analytical re- 
sults. 

The validity of Eq. Q and of our method to approx- 
imate Pxim) at low temperature, have been tested for 
a wide range of values of a (even for short range inter- 
actions) and a detailed study will be reported elsewhere 



VI. CONCLUSIONS 

In conclusion we propose a general method to deter- 
mine the anisotropic energy barrier in spin systems. The 
barrier, which determines magnetic decay times, can be 
computed from the topological non-connectivity energy 
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threshold in the corresponding isolated systems. Our 
analysis shows that for long range interaction two dif- 
ferent regimes can be identifyied, coherent rotation and 
nucleation. In the first regime we have predicted and nu- 
merically confirmed that magnetic decay time depends 
exponentially on the volume of the particle. Neverthe- 
less, in the regime of nucleation, the magnetic barrier de- 
pends on a suitable power of the volume, v 2 ~ a ^ d , where 
a < d is the range of interaction and d is the lattice di- 
mension. These facts have two remarkable consequences. 
The first one is that magnetic decay times depend ex- 
ponentially on the volume of the particle and not on its 
cross sectional happens for short range interac- 

tions. The second one is that, the fast growth of the 



magnetic barrier with the volume, gives rise to the pos- 
sibility to observe long-lived metastable ferromagnetism 
in finite systems at high temperature which is ruled out 
for short range interaction. 
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